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Abstract 

We implement relativistic BCS superconductivity in = 1 supersymmetric field 
theories with a U (1)_r symmetry. The simplest model contains two chiral superfields 
with a Kahler potential modified by quartic terms. We study the phase diagram 
of the gap as a function of the temperature and the specific heat. The supercon- 
ducting phase transition turns out to be first order, due to the scalar contribution 
to the one-loop potential. By virtue of supersymmetry, the critical curves depend 
logarithmically with the UV cutoff, rather than quadratically as in standard BCS 
theory. We comment on the difficulties in having fermion condensates when the 
chemical potential is instead coupled to a baryonic U{1)b current. We also dis- 
cuss supersymmetric models of BCS with canonical Kahler potential constructed 
by "integrating-in" chiral superfields. 
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1 Introduction 

Superconductivity is a common phenomenon that arises whenever there is spontaneous 
symmetry breaking (SSB) of a local U{1) symmetry. BCS is a particular theory realizing 
SSB, where one starts with a theory with a local U{1) symmetry and quantum effects at 
finite chemical potential generate a SSB vacuum by fermion condensation. The IR choice 
of vacuum can be described in terms of an effective Landau- Ginzburg theory (which can 
be derived from BCS). The low energy excitation spectrum can be described in terms of 
a Landau liquid, where the excitations are fermions. 

The plan of this work is to investigate the extent to which the dynamics of relativis- 
tic BCS theory [HE], with its usual features, can be implemented within the context of 
jV = 1 supersymmetric field theory. A supersymmetric model for chiral symmetry break- 
ing produced by fermion condensation at zero temperature and zero chemical potential 
was discussed in |3]. This model was generalized in [4j to incorporate BCS type supercon- 
ductivity, but the construction uses explicit supersymmetry breaking terms -therefore the 
Lagrangian does not describe a supersymmetric theory (in addition, it involves approxi- 
mations where some terms of the Lagrangian need to be neglected). To our knowledge, 
there has been no discussion in the literature implementing BCS superconductivity in 
supersymmetric theories. 

BCS requires the introduction of chemical potential for the fermions and in super- 
symmetric theories this leads to some obvious problems. Consistency demands that this 
chemical potential be coupled to a (non-anomalous) U{1) current. For a baryonic U{1)b 
symmetry, in supersymmetric theories, this can only be done in a consistent way by si- 
multaneously introducing the same chemical potential for the scalars. But charged scalar 
fields in the presence of chemical potential can run into problems of Bose-Einstein (BE) 
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condensation when the chemical potential becomes greater than the mass. The one-loop 
potential becomes ill-defined due to divergences. Adding a mass term to the superpo- 
tential does not circumvent this problem because the requirement of existence of Fermi 
surfaces, due to the relations between mass parameters implied by supersymmetry, is al- 
ways correlated to the appearance of BE condensation. We will evade this problem by 
coupling the chemical potential to a f/ (1)^ current and considering models where the light 
scalars have vanishing f/(l)_R charge. 

Another approach (used in OE]) is to compute thermodynamics quantities for the 
theory on x S^, where the three-sphere has radius R and the scalar fields have a mass 
equal to 1/R. Then the free energy of the system can be computed in a certain regime 
of parameters, typically, for chemical potentials which are lower than (or equal to) 1/R. 
However, we will see that in this approach the scalar mass scale cannot be separated from 
the Fermi energy. Although this does not completely preclude the construction of models 
with fermion condensates, it nevertheless implies that any model of this sort will be on 
the verge of producing BE condensation by a slight modification of parameters. It also 
implies that the vacuum dynamics will be governed not only by fermions near the Fermi 
surface but it will also be strongly affected by the scalar field fiuctuations, which in some 
cases can be dominant. 

Some previous studies of phase transitions in supersymmetric field theories have not 
found superconducting phases. In particular, in the free energy for = 4 super 
Yang-Mills theory was computed in detail for two particular values of the chemical po- 
tential, = or fii = 1/R, i = 1,2,3 (associated with f/(l) x f/(l) x f/(l) c S0{6)r). 
However, so far no sign of a superconducting phase transition was found despite all 
rich phenomena that seem to be taking place on the gravity side at strong coupling [7]. 
Another detailed search for U{1) breaking transitions was carried out in [8] for J\f = 4 su- 
persymmetric Yang-Mills theory coupled to a single massive fundamental-representation 
M = 2 hypermultiplet, but no evidence of any instability was found. Since any theory 
with higher supersymmetry can be viewed as a particular M = 1 supersymmetric field 
theory, it is convenient to use the M = 1 framework to provide a general picture of the 
conditions under which U{1) SSB can arise by BCS fermion condensation. 

This paper is organized as follows. In section 2 we briefiy review relativistic BCS 
theory, as this will provide the basis for the construction of a BCS theory in the su- 
persymmetric case. In section 3 we construct M = 1 supersymmetric Lagrangians with 
quartic fermion interactions. In section 3.1 we discuss a model with Fermi surfaces but 
with problems of BE condensation. We show that the one-loop potential is ill-defined 
even in 5*^ x 5''^ as soon as Fermi surfaces appear. In section 3.2 we present a simple 
example of an = 1 supersymmetric theory that exhibits BCS superconductivity with 
no problems of BE condensation arising from the scalar sector. In section 4 we consider 
models with canonical Kahler potential, and discuss the difficulties in implementing BCS 
superconductivity in these type of theories if all fields are dynamical and interactions ex- 
clusively arise from the superpotential. Section 5 summarizes the results and Appendices 
A and B contain further details of the calculations. 
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2 Relativistic BCS theory 



It is useful to briefly review the main features of relativistic BCS theory. Here we will 
work with global U{1) symmetries, so in this sense we will be discussing superfluidity, 
although transport properties are similar in both cases. In relativistic BCS theory, one 
has the effective Lagrangian [2] 

2 

C = '-{tp-f^d^iJ - d^^^^ij) - mijij + fiij^ij + ^(^cTs^)^^^) • (2-1) 

The U{1) symmetry ensures fermion number conservation, which allows the introduction 
of the chemical potential in the usual way. The Lagrangian is not renormalizable, the 
four-fermion interaction typically represents an irrelevant operator, but the dynamics of 
BCS superconductivity is such that for fermions which are close to the Fermi surface this 
attractive, four-fermion interaction becomes strong. At weak coupling the scaling dimen- 
sion of the fermionic fields must be very close to that of the 3/2 for a free field. Hence, on 
dimensional grounds, the interaction term is irrelevant in the IR. Naively it would seem 
that this theory cannot lead to any interesting IR physics. The phenomenon that actually 
takes place is explained in [9]. The key observation is that, in the presence of a chemical 
potential there is a Fermi surface which can change the naive scaling dimensions for the 
operators in such a way that the otherwise irrelevant interaction becomes indeed marginal. 
This is the seed for the possibility of a non-trivial IR physics such as superconductivity. 

Technically, to study this system, one considers the Euclidean theory at finite temper- 
ature and performs a Hubbard-Stratonovich transformation, one introduces the auxiliary 
field A{x) and the Lagrangian is then replaced by 

1 i - ■ 

+ ^1^1' - ^ [At (^,75^) + A(^,75^)t] . (2.2) 

The Lagrangian becomes quadratic in the fermions, which can now be integrated out 
explicitly giving rise to an effective potential for A. The fermion energy eigenvalues are 

u± = ^/Mp}±fiy + \A\^ , ^^0 = V^T^ , (2.3) 

where ± stands for particles and antiparticles. The one-loop effective potential Q is then 
obtained as usual by adding to the classical potential, ^|Ap, the thermal contribution 

-| /^(log(l + e-/^--(-)) + log(l + e-/'--^(^))) , 
plus a (Coleman- Weinberg) contribution that survives at zero temperature, 

/ (0j(2^o(p)-a;-(p)-c^+(p)) ■ (2.4) 
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is the thermodynamic potential of the grand canonical ensemble. The integral over 
momentum for this contribution is divergent. Since the theory is not renormalizable, one 
must restrict to energies below a cutoff A ("Debye" energy). The cutoff as usual represents 
the energy where new physics emerges. At low temperatures, the dominant contribution 
then arises from frequencies ujq near /i and the contribution of the antiparticle -represented 
by the terms with a;+- can be neglected (we assume > 0; if /i < it is the antiparticle 
contribution the dominant one). Since > 0, the existence of a Fermi surface at a finite 
momentum pp requires ji > m, where Pf is defined by the condition ^/p]? + rn? = j2. As 
a result the system has a Fermi energy represented by the chemical potential /x. If A = 0, 
this represents the Fermi energy in the usual sense, at zero temperature fermions would 
occupy energy levels with uJo{p) < /i. However, in this system, A is spontaneously turned 
on below some critical temperature. When A is not equal to zero, there is a fermion 
condensate and the energy eigenvalues u}±{p) do not vanish at any value of momentum. 
At low temperatures, the dominant contributions come from the region where uj-{p) has 
a minimum value. The location in momentum space of this minimum value defines the 
concept of Fermi surface in more general situations. For this system, this still occurs at 
oJo{p) — /i, though we will see that in more general systems the Fermi surface location 
can be shifted when A is turned on. 

The instability leading to A 7^ and thus to the formation of the fermion condensate 
appears when the coefficient of the (9(A^) term in the complete expression for the one- 
loop effective potential changes sign. Let us examine the conditions under which the 
fermion condensate appears. Expanding the full one-loop effective potential Q including 
the Coleman- Weinberg and thermal part in powers of A, one finds 

01 -^Vl f\ .f tanhilPMp)-^)) tanh(|^(a;o(p)+/.)) x\ 

The equation ^^|^2 = determines the critical temperature for the formation of the 
fermion condensate. As long as > m this equation always defines a finite critical 
temperature for any value of g, m, /i and A ^ m. The gap equation A = A(T') is 
obtained by differentiating the one-loop effective potential with respect to A. We find 

1 ^ ^ '^^-<':-^" + i'^^^-f'^n . (2.6) 

2ir'J„ '^'^y u)_(p,A) t^+(p,A) / ^ ' 

This gives the usual critical curve for a second-order phase transition for the order pa- 
rameter A as a function of the temperature. There are standard approximations that 
one can do. The second term inside the integrand comes from the antiparticle and can 
be neglected as ou^ <^ uj+ near the Fermi surface. In doing so one connects with the 
expressions of the non-relativistic case. Near the Fermi surface one may also approximate 
the factor p^ in the numerator by /j,^ — w? . 
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3 Supersymmetric BCS 



Let us now try to design a supersymmetric Lagrangian which incorporates these basic 
features. We are interested in a supersymmetric theory with a global baryonic U{1)b 
or U{1)r symmetry which undergoes spontaneous symmetry breaking. In general, U{1) 
symmetry breaking is easy to achieve by a suitable choice of the superpotential W. How- 
ever, here we are looking for a BCS type mechanism, where the breaking is caused by 
fermion condensation triggered by quantum effects. Af = 1 supersymmetric models with 
a canonical Kahler potential do not contain any quartic fermion interaction for any choice 
of superpotential W (see section 4). Quartic fermion interactions arise by means of the 
following choice of Kahler potential: 



We would like to construct a supersymmetric BCS theory with Dirac fermions, and 
in A/" = 1 supersymmetric theories this requires at least two chiral superfields (a single 
chiral superfield describes a Weyl fermion). The simplest theory consists on two chiral 
superfields X and Y with the Kahler potential 



The coupling g could in principle be different for the interaction terms involving X and 
Y superfields. One could also add, for example, a term X'^XY'^Y (used in [3]). However 
we shall consider the above simple choice which already illustrates the essential points. 

3.1 Chemical potential for U{1)b 

We first consider the J\f = 1 supersymmetric model defined in terms of two chiral super- 
fields with Kahler potential (13.21) and superpotential: 



This gives masses to scalars and fermions. It will be shown that this model is not suitable 
to implement BCS mechanism in supersymmetric theories. The model will illustrate the 
typical problems that one has to deal with. 

We first consider the Lorentzian theory on R^. For the finite temperature theory, 
we shall later consider the Euclidean theory on x R^, and eventually on x S^. In 
components, the Lagrangian reads 



is:($,$t) = $t$ + /($t$)2 _ 



(3.1) 



K{X, F, X\ Y^) = X^X + Y^Y + g\X^Xf + g\Y^Yf . 



(3.2) 



W = mXY . 



(3.3) 




^^(^.^,.)(^t^t) 

1 + V|0^.|2 



(3.4) 



(3.5) 
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Here we use rj^i, = diag(l, —1, —1, —1) and as usual a^'"" = (1, — cr*), where a* are the 
Pauh matrices (we will follow the notation of [TP]). 

Note the presence of the (non-renormalizable) quartic fermion interaction. The cou- 
pling constant g has dimension of length. The choice of sign of g"^ was made in order to 
have the same type of interaction as in BCS. We have checked that the opposite sign does 
not lead to fermion condensation by quantum effects. For g"^ < there is no consistent 
solution to the gap equation for the vacuum condensate. The effective potential is un- 
stable and cannot be consistently minimized in the one-loop approximation. Therefore in 
what follows we assume g^ > 0. 

We need to introduce a chemical potential and consistency demands that this is cou- 
pled to a conserved non- anomalous U{1) current. The superfields X and Y carry opposite 
f/(l) charge so the baryonic U{1)b current is non-anomalous. Turning on a chemical po- 
tential corresponds to turning on a background U{1)b gauge field component Aq = /i. 
In order to have a Lagrangian quadratic in fermion fields, one can make a Hubbard- 
Stratonovich transformation in the component Lagrangian by introducing two auxiliary 
fields, A^, Ay, 

£5 = (1 + - ^''tu 12 - 9'{i + vi0.niA.p +ix^y) m 

l + Ag^\(f)x\^ 

+ ( {y^J^^ + 9'^^) (^.^x) - Irni^^^y + h.c)j + {x ^ y) , (3.7) 



where = d„ — iq^S^o (with no loss of generality one can set the X U{1) charge qx = 1, 
as it can be absorbed into a redefinition of /i; in this way Y has charge qy = —1). The 
Lagrangian has now become quadratic in the fermion fields, no quartic fermion interaction 
is left. As a result, the functional integral over fermions can be directly performed. 

Next, we expand the scalar fields around their VEV's, (p = v + (p, and retain only up 
to quadratic terms in the scalar fields (we assume real v). We find 

Cs = {l + Ag\l)d,^ld>^^^ + Ag\l - .""^'ff', ) ivl + vf) (3-8) 



+ n , A 9 o^2 if^'fy + fxfy + VxVy + VxVy 

(1 + Ag^vl) 

+ ((1 + 16.^.)/^ -4. IA.I ^^^-^-^ jl^.l 

+ i\x{\ + 'ig^vW^ldt^x - VxdtVl) - Aifig\l{ipldt(pl - <Pccdt<Px) + {x ^ y, fi ^ -fi) 



(3.9) 

U = t{l + Ag^lMa^d.^x) + /i(l + Sg'vlMa'tlJx) 
' 2mg'^VxVy 

l + Ag^vl ' " ~V 2" 



+ ( ( .,.22 + 9^^x ii^xi^x) - -^i^xi^y + h.c. ) + (x ^ y, /X ^ -/i) (3.10) 
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with 



2 2 

m vz 



^ci = YTlghP + + A.|2 - + (x o y) . (3.11) 

To have canonically normahzed kinetic terms, one can redefine fields as follows: 

a/1 + Ag^v^ a/I + 4g'2^;2 

Integrating over 99, Lp* leads to a one-loop potential depending on v, A, n, m. Since 

the model is not renormalizable (just like BCS) integrals will be regularized by a momen- 
tum cutoff, representing a "Debye" energy where new microscopic physics appears. 

We proceed as follows. CaUing Os and to the resulting 4x4 scalar and fermion 
matrices for the quadratic terms in momentum space, we shall write the determinants as: 

4 4 

det O5 = n - , det Of = n " ' (^.13) 

1=1 1=1 

where 

ujsi= <^SiilJ',\p\,9,m,V:c,Vy,A^,Ay) , ujpi = (^FiilJ',\p\, g,m,v^,Vy, A^, Ay) . (3.14) 

The expressions for Os and Op are shown in appendix B. The eigenvalues for the fre- 
quencies have complicated expressions when and Vy are non-vanishing. The strategy 
is to look for non-trivial minima at = Vy = with A^., A^ 7^ 0, assuming them to be 
real. Next, we shall check that the one-loop effective potential is locally stable in Vx and 
Vy directions, a property that will be ensured by the presence of a mass term. 

When Vx = Vy = the scalar and fermion quadratic terms greatly simplify. At this 
point, we find the following eigenvalues for the frequency. 

ous 1,2 = V^9^^l + m'^+p'^ ± /X , 

3,4 



^4^4 A2 + m2 + p2 ± ^ ^ (3.15) 



ul 3,4 = 2g^Al + 2g^Al + fj,^ + + ± S^, (3.16) 

£:± = 2^/.2 (m2 +p2) + ^8 - A2) 2 + ^4 (^2 + A,)' ± 2/.p (A^ - A^)) (3.17) 
For configurations with Aj; = A^ = A, the fermion frequencies become 



\ 



Ip'^ + + 4g^A^—±^\ + 4g^A^ [l - — ) . (3.18) 



On the other hand, for A^; = 



— Aj, = A, we find 



Up = \/{Vp^ + ± /^)^ + 4^^A2 . (3.19) 

This is the same dispersion relation as in the relativistic BCS system of section 2. This 
might suggest that BCS mechanism can be implemented in a similar way. But the presence 
of charged scalars demands some care. We first need to identify the Fermi surfaces. For 
A^. = Ay = 0, they lie on the region where up 2,4 vanish, i.e. at 

^Jpl + m2 = fi . (3.20) 

As in the standard relativistic BCS case, the existence of a Fermi surface would require 
fi > m. However, in the present supersymmetric system we cannot set fi > m because 
the scalar contribution to the thermal partition function 

^E/ >°^(l (3.21) 

is ill-defined, because us 2,4 become negative below some momentum. The system presents 
BE condensation, the occupation number of scalars with zero momentum goes to infinity 
as /i approaches m from below. This spoils the BCS mechanism. 

One possible approach to elude this problem while maintaining supersymmetry is to 
put the theory on 5"^ x 5"^. Because the scalar field couples to the curvature (see e.g. [TT]). 
this will provide an extra mass term for the scalar fields, which might allow for regions 
in parameter space with Fermi surfaces, and without problems of BE condensation. The 
mass term, when the i?-charge of the scalars is one, is now of the form 

[m' + R-') (0:0, + <j);<Py) , 

where R is the radius of the three-sphere. The scalar contribution would be negligible if 
one could assume that 1/R> A. However, having put the theory on S^, the integral over 
momentum is replaced by a discrete sum originating from the Kaluza-Klein modes of S^. 
This replacement is achieved by 

Scalars : / — > 1(1 + 2)R-^ 

Fermions : — > {l + l/2fR-'^ (3.22) 

with / = 0, 1,2, . . .. One must also take into account the degeneracy: for scalars, df = 
{I + 1)^; for fermions, df = l{l + 1). In particular, for the fermions, / = does not 
contribute. For the scalars, in addition we must add the mass term R^^. This is 
effectively incorporated by the replacement 

/ /(/ + 2)R'^ + R-^ = {l + ifR-^ (3.23) 
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These formulas show that one cannot assume 1/R > A, since such cutoff would leave no 
excitation in the system. Therefore it is not possible to separate the scalar mass scale 
from the Fermi surface. 

In order to see if the system can have Fermi surfaces, we need the detailed form of the 
Lagrangian on S^. This depends on the R charges of the fields. We denote by q the R 
charge of (px so that the charge of 0^ is 2 — q. From the expressions given in [11] , we find 



- 2) , „ g - 1 2^ I , 12 , , 2 - 1) , 2g 



m^\(f)y\ 



2u 12 + Vl0xr)|A.|^ + (x^2/,/i->-/i,g-^2-g) (3.24) 

1 + 45(^10.1^ 

+ + /i) (^i^V.) + V + 2/i) |0.P(^icrV.) 

-(ifi^ifr-^S^^^^^^^^ 

- -miijjx'iljy + tpli/jl) + {x ^ y,iJ,-)- -iJ,,q 2- q) (3.25) 

We shall demand that in the unbroken phase the theory has well-defined thermody- 
namical potentials. So we begin by considering the case A. = A^, = 0, Vx = Vy = 0. We 
will now see that BE condensation is inevitable in this case, which is sufficient to rule out 
the model. Consider first the case q = 1, i.e. the U{1) charges of X and Y are equal to 
1. The scalar contribution is now given in terms of the frequencies 



Us 



V(/ + 1)2/^-2 + m2 ± /i , / = 0,1,2,... (3.26) 



If both X and Y had the same baryon charge, the Fermi surface would just be determined 
by the replacement fl3.22p in the fiat expression f l3.20p . and shifting the chemical potential 
by /i — )■ yU + l/(2i?). As X and Y have opposite baryon charges, this is more involved. 
By explicitly computing up from the above Lagrangian, we obtain that the Fermi surface 
cji? = is at 

yJl^R-"^ + m? = ^i, If = 1,2,... (3.27) 

For a given choice of Ip, one can determine Substituting in the lowest {I = 0) scalar 
frequency, we see that the scalar frequency cannot be positive as long as = 1, 2, . . ., 



VR-^ + m2 - ^/2,i?-2 + ^2 < . (3.28) 

Therefore, even on S^, it is not possible to separate the Fermi surface from the region of 
BE condensation. The underlying reason being that the extra mass term for the scalar 
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provided by the coupling to the curvature of the space is of the same order as the quantized 
fermion momentum values. The same problem arises for any choice of q. 



3.2 A simple supersymmetric BCS model: 
Chemical potential for U{1)r 

Let us now consider an A/" = 1 supersymmetric model with two chiral superfields X and 
Y with Kahler potential given by (13. 2p and superpotential W = Q. The Lagrangian 
has a U{1)r symmetry for arbitrary U{1)r charges of the X and Y superfields. It is 
convenient to consider the U{1)r symmetry under which scalars (px and (f)y are neutral, 
so that fermions ipx and ipy have charge —1. The advantage of this choice is that we can 
avoid problems of BE condensation even in R^. Note that with this charge assignation 
the U{1)r symmetry is anomalous. However, this can be easily cured by adding to the 
theory free superfields with canonical Kahler potential with the required U{1)r charges 
to cancel the anomaly. For example, one may add Zj, i = 1,2 with i?-charges R{Zi) = 2 
so that tpzi, '4^Z2 have charges +1. The scalars in Zi would then couple to the chemical 
potential and may undergo Bose-Einstein condensation. However, this sector is completely 
decoupled and therefore does not participate in the thermodynamics governing the X, Y 
sector. 

The component Lagrangian with chemical potential included can be obtained from 
the previous case, f l3.6p . f l3.7p . by setting m = 0, vanishing U{1) charges for the scalar 
fields (which amounts to replace covariant derivatives of the scalar fields by ordinary 
derivatives) and taking into account that fermions ipx and ipy now have the same charge 
— 1. The quadratic Lagrangian for the fluctuations (after expanding around expectation 
values) is given by 

Cs = d.^ld^^x + d.^ld^^y - Y^^W.? - n^l^^l' ' ^3.29) 

where we have rescaled the fields to have canonical kinetic terms. The classical potential 
is given by 

= {^g'vl + 1) \Axf + ix^y). (3.31) 
The equations of motion for A^,, Aj^ give (setting the scalar fluctuations (px,'^y -> 0) 

A _ ^i^i A _ ^^^^ 



[1 + Ag^vl)^ ' (1 + Ag\, 



Ax, Ay have both U{1)r charges equal to 2. Vacuum expectation values for them thus 
spontaneously break U{1)r and represent a measure of the fermion condensate. 
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By proceeding in a similar way as in the previous case, we now find the following 
frequencies for scalars and fermions 



S 1,2 



2^ VA 

p + 



F 1,2 



1 + Ag^v. 



2„,2 ' 



S 3,4 



2^ 

J9 + : 



%^A2 



F 3,4 



%^A2 



(1+Vt;2)' 



(3.33) 
(3.34) 



Here we have chosen real A^., Ay, as one-loop potential depends only on their moduli. We 
stress that these simple dispersion relations are a consequence of the extreme simplicity of 
this supersymmetric model; generic models (even with simple superpotentials) typically 
lead to very complicated eigenvalues for the frequencies. 

In the present case, the dynamics of the X and Y fields are decoupled. It is clear 
that the same configuration that minimizes the one-loop potential in the x direction also 
minimizes the one-loop potential in the y direction. Therefore with no loss of generality 
we set Vx = Vy = V and A^. = A.^ = A. 

The complete one-loop thermodynamic potential is given by 



n = 2g' (1 + Ag\^) A^ + -L ^ dp Ulog 



sinh 



/3 



Ip2 + 



Ag^A^ 
1 + Ag^v^ 



log 



cosh^^/(p + /i)2 



Ag^v 



2,,2-\2 



-log 


cosh^^ 







4C/4A2 



4^2^2)2 

(3.35) 



When the vacuum lies at A 7^ 0, then f = is a local minimum. When the vacuum lies 
at A = 0, then there is a flat direction in v, because in this case the frequencies do not 
depend on v. This is confirmed by the evaluation of the one-loop potential. 
The gap equation A = A(T) is determined by the equation 



dA 







(3.36) 



This gives, when v = 0, 



dp p^ 



27r2 



tanh (^i/3v/VA2 + (p - tanh Ag^A^ + {p + ^if 



^Ag^A^ + (p - /i)^ 



^Ag^A^ + (p + iif 



2coth i/3v/4^W 



p" 



^Ag^A^ + p2 



(3.37) 



The gap equation can be compared with the gap equation (12. 6p of the relativistic BCS 
system of section 2. One difference is that now scalars and fermions have zero mass, since 
a mass term mXY would not be consistent with scalars neutral under f/(l)^. The second 
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and more fundamental difference is given by the scalar contribution -represented by the 
second line in the above equation- that we analyze in what follows. 

An important consequence of supersymmetry is that the critical curve A(T) now 
depends logarithmically with the cutoff: for large p, the integral in fl3.37p behaves as 
follows: 2 A 2 

9 f dp log A. (3.38) 

TC^ J p 

If the scalar contribution is removed, like in non-supersymmetric BCS, one has instead 

dpp^ . (3.39) 



,2 M 



7r2 

Obviously, a logarithmic dependence with the UV cutoff is a desirable feature, since the 
thermodynamics becomes much less sensitive to the underlying microscopic physics. 

At the same time, the IR physics produced by the scalar sector has a striking effect: 
the superconducting transition becomes first-order, instead of second-order, as it would 
be in standard BCS. The IR physics of the scalar sector is important at the onset of the 
transition, where A is small. To see the nature of the transition, we need to compute 
d^/dT. This can be obtained by differentiating the gap equation dfl/dA with respect to 
T. Writing the gap equation in the form 1 = /(A^,T), one has 

dA 1 ^ 

dT 1^ Jf- ^ ' 

In a second-order phase transition, dA/dT is singular at the critical temperature, where 
A = 0. This is because df /dT and df/d{A'^) are regular at A = 0. While the scalar 
contribution to the one-loop potential is regular at A = 0, its second derivative with 
respect to A^ has a singularity near A = originating from the region near p = 0. We 
have 

^/ (3.41) 



a(A2) 7r2 7o (4AV+P^)^ TT A 

As a result, dA/dT is now finite at A = 0. The superconducting phase transition is 
therefore first-order. This significant change coming from the p = region would obviously 
not take place if the scalar field was massive. In such a case, the phase transition would 
still be second-order. But, as explained, in the present model it is not possible to add a 
mass term. 

A = A(T) is shown in fig. [T]for different values of the chemical potential. We see that, 
as the chemical potential gets smaller, the transition approaches to a second-order phase 
transition. In general, the scalar field has the effect of decreasing the critical temperature 
with respect to the relativistic BCS case. 

In an interval of temperature, T^i < T < Tc2, there are three branches -characteristic 
of first-order phase transitions- corresponding to three solutions of the gap equation: the 
trivial minimum at A = 0, a maximum and another minimum at higher A. These are 
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exhibited in fig. |2l showing the one-loop potential at different temperatures. We see how 
the non-trivial maximum and minimum are created as the temperature is lowered below a 
certain critical value (Tc2 ~ 4.2 in fig. [21 see also fig. [T]with /i = 10). At a temperature 
Tc, the non-trivial minimum becomes degenerate with the minimum at A = (in fig. |2l 
this occurs at Tc ~ 4.09). In the interval Td < T < Tc, the symmetric vacuum A = is 
metastable. Below ~ 3.55, the symmetric vacuum A = becomes unstable and the 
only minimum of the potential is the SSB vacuum at A 7^ 0. 




1 2 3 4 5 0.0 0.2 0.4 06 0.8 1.0 1.2 14 0.00 0.05 ttlO 015 

T T T 



Figure 1: Gap A vs. T at = 10, 4.5, 3 = 0.5, A = 40). 




Figure 2: Potential showing generation of minimum at A 7^ 0. Different curves (from top 
to bottom) at T = 4.2, 4.17, 4.09, 4, 3.55 (with g = 0.5, /i = 10, A = 40). 

The gap at T — )■ can be determined analytically. In this limit, the hj^erbolic tanh 
and coth become 1 and the integrals can be easily computed. Assuming A ^ gA, fi, we 
find the following result 

l-^(2>ogf^V3l. (3.42) 



27r2 V \g'^A{0) 

i.e. 



TT 3 

A e~^~2 

A(0)^ . (3.43) 
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It should be noted that our normahzation for A is different from the normahzation of A 
of section 2. One can easily go back to the normalization of section 2 by the substitution 
2(7^ A — 7- AsusyBcs- Taking into account this, we are finding 



3 



|A(0)susyBcs| ~2Ae ^ . (3.44) 

This is similar to the standard formula in BCS theory with the identification /i^ = 7r^A^(0), 
where iV(0) represents the electronic density of states at the Fermi energy. It is useful to 
reproduce the non-supersymmetric formula in the present context. We consider the gap 
equation fl3.37p . dropping the scalar and the antiparticle contribution. As explained in 
section 2, the latter is negligible near the Fermi surface. We obtain the equation 

o2 /-A I 
1=-^/ dpp'^ , (3.45) 

This formula diverges quadratically. It does not give a sensible result if the 'Debye energy' 
A is far from the Fermi energy. The standard procedure involves an ad hoc approximation, 
where in the numerator is replaced by fi"^ and the integral is done in the interval 
\p — n\ < A. This leads to the result 



|A(0)bcs| ~2A 3^ . (3.46) 

1 — e s^f^ 

In the supersymmetric case, because quadratic divergences are canceled between fermion 
and scalar contributions, the analog result (13.44^ follows by performing the integrals in 
the gap equation exactly, with no need of the ad hoc approximation — > /i^. 

Let us now study the specific heat. It is instructive to examine the different contribu- 
tions to the thermodynamic potential closely. Consider first the symmetric phase A = 0. 
As we shall be interested in derivatives with respect to the temperature, we can subtract 
the Coleman- Weinberg contribution, so that integrals are convergent. We write 

r2 ^scalar ~l~ ^electron ~l~ ^positron ; (3.47) 

where (after integration by parts) 

2 e~T 

1 r°° e-^ 

Jo 1 + e T 

^positron|^^Q = ~7r^ / dp p^ ZeIr ' (3.50) 

Jo 1 + e T 
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We get 

2 



as expected, since, when A = 0, the scalar contribution describes a relativistic boson 
particle (there is an extra factor of 4 as compared with the usual single scalar contribu- 
tion, because we have two complex scalar fields ipxi (fy). The integrals for the fermion 
contributions can be computed analytically in terms of polylogarithmic functions. 

Let us now compute the difi^erent contributions to the specific heat. The entropy and 
specific heat are given by the familiar formulas 

Note that specific heat here is defined as a partial derivative at constant /i, instead of 
constant charge density p. In the present case, we find this to be a more sensible quantity, 
since the scalars are neutral and would not contribute to the constraint p = dfl/dp,. 
We obtain 

Cscalar\^^Q = ) (3.53) 

as usual for relativistic bosons. Consider now the fermion contributions. At large T, the 
dependence on p, disappears and one gets the usual behavior of a relativistic fermion 



^electron | q ^positron | 



= , for T » . (3.54) 



At low temperatures, CposUron is exponentially suppressed, CposUron ~ e"'^^'^ . For the 
electron, the integral picks the main contribution near the Fermi surface, p p., and one 
gets the usual linear behavior for the electronic specific heat at low temperatures 

C = Ceiectron] ~ ■ (3.55) 

Let us now compute the full c{T) including the region T < Tc2 where A 7^ 0. We use 
the notation e — A"^. Then 

<T) = -T i—^ + ) , ^ = - w , (3.56) 

where we used the fact that A(T) is defined by dfl/dA = 0. The resulting c(T) is shown 
in fig. 3. At T < Tc2, it exhibits the expected exponential suppression due to the gap. 
At the critical temperature, we notice the characteristic discontinuity of first-order phase 
transitions. In first-order phase transitions the entropy may experience a finite jump at 
the transition, leading to an infinite jump in the specific heat. For Tc2 < T < 0{pi), we see 
the linear behavior coming from the electron contribution. Finally, at high temperatures, 
it exhibits the behavior shown above in terms of analytic formulas. 
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12 3 4 5 6 7 

Figure 3: Specific heat as a function of the temperature (with g = 0.5, /i = 10, A = 40). 

4 Models with canonical Kahler potential 

Let us now consider the possibility of implementing a BCS mechanism in models with 
canonical Kahler potential. Fermion condensation now needs to be triggered by inter- 
actions contained in the superpotential W. A way to find an appropriate model is by 
starting with the previous model with non-canonical Kahler potential and "integrating 
in" some fields. One simple example is the model defined by the superpotential 

K = X^X + Y^Y, W = moZ{X - gY^) . (4.1) 

The chiral superfield Z has no kinetic term and can be integrated out. The equation 
for Z sets X = gY"^ . Substituting this relation into the Kahler potential, one finds a 
new Lagrangian with K = Y'^Y + g'^iY'^YY , which contains quartic fermion interactions, 
as desired. Thus this model is equivalent to the model considered in section 3 (after the 
addition of another superfield similar to Y). The model f l4.2p represents a supersymmetric 
analog of Hubbard-Stratonovich Lagrangian, with the chiral superfield X playing the role 
of A. Note that there is a U{1)b as well as a f/(l)_R symmetry. One can choose a U{1)r 
charge assignment such that X and Y are neutral and Z has charge equal to 2. 

Another example is given by a supersymmetric model containing three chiral super- 
fields X, y, Z, with potentials 

K = X^X + Y^Y , W = moZ{X - gY'^) + MZ'^ . (4.2) 

Once again, Z has no kinetic term and can be exactly integrated out. We end up with 
the effective model 

K = X^X + Y^Y , W = m(X - gY^f , m = -—^ . (4.3) 

It should be noted that this model is not equivalent to (14. ip . If one ignored the kinetic 
term for X, the equation for X would set X = gY"^ (the same result of course follows by 
working in terms of component fields). But it is not legitimate to ignore the kinetic term 
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X't^X as this is precisely the term that will generate the effective four-fermion interactions 
and the dynamics that we wish to study. Therefore, understanding the vacuum dynamics 
requires a detailed examination of the model. 

In what follows we shall investigate the three models with superpotential (14.11) . (14.21) 
and (14.31) . but now with canonical Kahler potential for all fields X, Y, Z. The models 
(14. 2 p and (14.31) have only f/(l)_R symmetry, under which Z and X have charge 1 and 
Y has charge 1/2. In all cases, the U{1) symmetries are anomalous, but this can be 
cured as before by adding a suitable set of free superfields. We would like to see if these 
models can reproduce a BCS mechanism leading to a fermion condensate, similar to the 
model of section 3.2, even in the case when all fields are dynamical. The basic Wess- 
Zumino component Lagrangian for the Af = 1 supersymmetric model with canonical 
Kahler potential reads 

C = d,<P*d^<P, + t4a^d,i^, - i {{d,d,W)^i^ii^, + {did,W%4i^^^i - \{d,W)A' ■ (4.4) 
In particular, the component Lagrangian for (14. 3 p is 

-4m2(l + 4/|0,p)|0.-^?0j|2 

+ {-mip^ip^ + 2mg{(f),j, - 3g(f)l)i!yi!y + Amg(j)yilj^iljy + h.c) . (4.5) 

Consider now the equation of motion for 0*. For constant fields, one would have 

2m(l + Ag^\ct>y?){ct>. - gct>l) = g^^ . (4.6) 

This equation connects the fermion bilinear V'y'^y with the expectation value of (px, in 
much the same way in the BCS Lagrangian the fermion bilinear is connected with A. In 
addition to ipy, its scalar superpartner (j)y also appears in the equation. If the term X'^X 
were removed from the Kahler potential, then the equation would just be (p^ = g(f)y. After 
the functional integral over ip, ip"^ and over the scalar field fluctuations is carried out, we 
are left with a one-loop potential depending on Vx = {4>x) and Vy = {4>y)- The expectation 
value Vx at f J, = is an order parameter and a measure of the fermion condensate. 

Standard BCS theory involves a Bogoliubov transformation, needed to diagonalize the 
Hamiltonian, which leads to the usual concept of quasiparticles. In the current system 
(14. 4 p with canonical Kahler potential, where (px ~ '^y'^yi Bogoliubov transformation 
is needed in virtue of the presence of Vxipyipy and v*ip^y'ip^ in the quadratic part of the 
Lagrangian. Diagonalizing the corresponding Hamiltonian involves a mixing of creation 
and annihilation operators. 

We summarize the results of the direct calculation for these models. 

The functional integral over fermions can be explicitly carried out as the Lagrangians 
are quadratic in fermions (the chiral superfield X represents the supersymmetric analog 
of the Hubbard- Stratonovich field A). The scalar fields are expanded around a vacuum 
value, e.g. (px = Vx + (Px, (py = Vy + ipy, (pz = Vz + (Pz- The conclusions are as follows: 
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• W = m{X — gV^y. The model has only U{1)r symmetry with {qx, Qy) = (Ij 1/2)- 
Chemical potential is introduced for this U{1). Like in the model of section 3.1 
with baryonic U{1)b chemical potential, scalar particles are charged. As a result, 
on one cannot have Fermi surfaces because they would overlap with the regions 
of BE condensation. Therefore, like in that model, we attempt to study the theory 
on a 3-sphere of radius R (relevant formulas are given in appendix A). At Vy = 0, 
there are Fermi surfaces at fiR = 2lp + 1, Ip = 1,2,.... For this value of /x, one 
finds a scalar frequency oJs = (1 + / — If)/R, I = 0, 1,2, .... This shows that it is 
not possible to have positive definite scalar frequencies for any I p. Thus the Fermi 
surface cannot be separated from the region of BE condensation where the thermal 
one-loop potential is ill-defined. 

• W = mZ{X — gY'^) + MZ^. This is a renormalizable model and can be viewed as 
an UV completion of the previous case. It has a U{1)r symmetry with uniquely 
determined U{1) charges, {qx,(lY,(lz) = (1, 1/2, 1), and a similar IR physics as the 
model W = m{X — gV^Y, with M playing the role of the UV cutoff A. 

• W = mZ{X — gV^). This model has a baryonic U{1)b as well as t/(l)i? symmetry. 
If chemical potential is introduced for U{1)b, again we find that the Fermi surfaces 
cannot be separated from the region of BE condensation (irrespective of the R-charge 
assignation). 

Consider now a chemical potential coupled to the f/(l)_R current. One can assign 
charges {qxAvAz) = (2 — 1 — g/2,g). Unlike the model of section 3.2, now 
it is not possible to have only neutral scalars. Nonetheless, for q = 2, i.e. when 
{Qx,(1y,<1z) = (0,0,2), in the unbroken phase = Vy = = 0, it is possible to 
have a Fermi surface without BE condensation even in fiat space. There is a Fermi 
surface at pp = fi. The scalar frequencies are 

LOs = {^/p'^ + m? ± 2/i, a/p^ + m?, \/ p"^ + m?, p, p} , (4.7) 

which are always positive definite for fi < m/\/2. The problem is that BE conden- 
sation reappears in an infinitesimal vicinity of = Vy = Vz = 0. After turning on 
Vz and Vx, the last frequency in f l4.7p becomes 

Us = \/ p'^ + Ag'^m?v1 — 2gm?Vx , 

which becomes complex at low momenta in the region Vx > 2gv1. Thus the model 
is not protected from BE condensation. 

It should be noted that there are many M = 1 supersymmetric models admitting 
superconducting phases where the favored classical vacuum is a SSB vacuum with t> 7^ 0. 
But this is a classical effect, it is not BCS superconductivity. In particular, it is not 
triggered by fermion interactions that become marginal near a Fermi surface. There are 
also models containing Fermi surfaces, i.e. regions in parameter space where the fermion 
frequencies vanish or take a minimum value at some fermion momentum but, for the 
reasons explained above, BCS type phase transitions do not occur generically. 
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5 Conclusions 



Understanding how superconducting transitions can take place in supersymmetric field 
theories in full detail is of great interest, in particular, to clarify how supersymmetric 
theories react upon the introduction of chemical potential. Other motivations include 
providing a field-thcorctical understanding of the possible mechanisms underlying holo- 
graphic superconductivity, and possible applications in real condensed matter systems 
containing fermion and scalar quasiparticle excitations. A chemical potential typically 
leads to Fermi surfaces for fermion fields, and to Bose- Einstein condensation for scalar 
fields. In a supersymmetric theory, fermion and scalars are combined with very specific 
couplings, which have a significant incidence in the radiative corrections that determine 
the one-loop effective potential. What is the impact of a supersymmetric combination 
of bosons and fermions on phase transitions in thermodynamical systems with chemical 
potential? What are the resulting phase diagrams? 

To address these questions, in this paper different models have been investigated in 
detail. The salient aspects of this investigation are as follows. The main obstacle to imple- 
ment BCS superconductivity in a neat way is, as expected, Bose-Einstein condensation. 
The model of section 3.1 exhibits the typical problems that arise. Introducing chemical 
potential for a baryonic U{1)b symmetry leads to the emergence of Fermi surfaces, but 
inevitably couples scalar fields to the chemical potential as well, since scalar fields have 
the same baryon charge as fermions. Supersymmetry prevents the scalar mass scale from 
getting separated from the fermion mass scale. Even if supersymmetry is spontaneously 
broken, the vanishing supertrace relation still implies that there must be light scalar fields. 
As a result, near the Fermi surface, the contributions to the thermal potential coming 
from scalar fields become ill-defined. In order to avoid Bose-Einstein condensation near 
the Fermi surface one can try to put the theory on of radius R. This gives an extra 
mass term 0{1/R) to the scalar field coming from the coupling to the curvature. How- 
ever, because the quantized fermion momentum is of order 1/R, the existence of Fermi 
surfaces requires chemical potentials that are also corrected by an amount of order 1/R. 
In this model we found that the Fermi surface cannot be separated from the regions of 
BE condensation. 

To circumvent these problems, in section 3.2 we proposed a specific supersymmetric 
model which reahzes BCS superconductivity in fiat space. It is based on a Kahler potential 
with quartic terms in the superfields. 

K = X^X + Y^Y + g^X^Xf + g^Y^Yf , (5.1) 

and W ~ 0. The chemical potential is introduced for a U{1)r symmetry under which 
the scalar components of the X and Y fields have vanishing charge. In this way Bose- 
Einstein condensation does not occur and the model can be studied directly on R''. This 
is presumably the simplest supersymmetric model for BCS superconductivity that one can 
construct, since it contains the minimum number of superfields to have a Dirac fermion, 
i.e. two chiral superfields, it can be studied on and the superpotential is 1^ = 0. 
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We found that the system has a superconducting phase transition below some critical 
temperature, produced by a fermion condensate. 

The equations determining the temperature dependence of the gap are very similar 
to BCS theory, with the main difference represented by the contribution coming from 
scalar fluctuations. One important effect of this contribution is a drastic reduction of the 
dependence on the UV cutoff from quadratic to logarithmic. Another effect due to the 
scalar superpartner is changing the character of the phase transition from second to first- 
order. As explained, the origin of the change is the contribution of low momentum scalar 
modes at small A. As the chemical potential is decreased, dA/dT at A = becomes 
large and the phase transition approaches to a second-order phase transition. 

In the zero temperature limit, the gap can be computed analytically. The resulting 
expression for A(0) is qualitatively similar to the standard BCS expression. The BCS 
formula is usually derived by integrating over momenta in a small neighborhood of the 
Fermi surface, otherwise one does not get a sensible result -for example, if the Debye 
energy was significantly larger than the Fermi energy, the formula would be different (for 
real materials, \ujc)ef,ye — w^l ^ up)- In the supersymmetric case, due to the cancellation 
of quadratic divergences, A(0) can be derived by integrating momenta in the whole range, 
from p = up to the cutoff energy, which may be much larger than the Fermi energy. In 
short, the supersymmetric expression turns out to be stable under variation of parameters. 
We have also computed the specific heat. In the superconducting phase, this exhibits the 
expected exponential suppression due to the gap. 

An interesting open problem is the construction of supersymmetric BCS models with 
local U{1) symmetry. In particular, this would permit the study of magnetic screening 
and more generally the response of the supersymmetric system to external sources. 
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A Thermal one-loop potential 

In the absence of chemical potential, the one-loop potential at finite temperature is com- 
puted by using the following formulas (see [12] for a recent discussion and references 
therein) : 

V = Vcl + Vcw + Khermal , (^-1) 
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where represents the classical contribution and 



1 M 

i 

^thermal = ^ (0 " ' ^^-2^ 



where i^^? = ^ Mf and the sum over i accounts for scalar {F = 0) and fermion 
{F = 1) degrees of freedom, and are the mass eigenvalues of the scalar and fermion 
mass matrices of the quadratic fluctuations, and A represents a UV cutoff. By expanding 
the logarithm in (1A.3P the integral over p can be carried out explicitly (see e.g. [12]). At 
sufficiently small temperatures, the contribution from massless modes dominate over the 
massive mode contributions, which become exponentially suppressed. This has the impor- 
tant implication that the SSB vacuum, which contains a massless mode -the Goldstone 
boson- will necessarily dominate the thermodynamics at low temperatures. 

When chemical potential is added into the system, the one-loop determinant is given 
in terms of eigenvalues of the frequencies. The complete one-loop potential is given by 
dAU) with 

i i 

where usi and upi are the scalar and fermion frequencies. We recall the origin of these 
expressions. In the Euclidean theory, time is periodic and one needs to impose periodic 
and antiperiodic boundary conditions for scalars and fermions, respectively. This leads 
to quantized frequencies, 

27rn 7r(2n-M) ^, 
^5 = , (^F = -^—^ , (A.6) 

with integer n. In particular, for scalars, one has to compute 

E ^ E + ^'s^) ■ (A.7) 

i n 

The sum over n can be easily performed by first differentiating with respect to Ug^. One 
finds 

^El°8((^)'+-«) = ilog.mh^ = f + i|og(l-e-''-) , (A.8) 
plus a constant which we discard. 
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A similar calculation for the fermions gives ^ + ^ log (l + e ^"^^^j . 

Consider now the construction of the one-loop potential for the theory on 5*^ x S^. As 
explained, once we put the theory on S^, the integral over momentum is replaced by a 
discrete sum over Kaluza-Klein modes of S^. The quantized momenta are as follows 

Scalars : — > l{l + 2)R-'^ , 

Fermions: / — > (/ + 1/2)^7?-% (A.9) 



with I = 0,1,2, .... Thus we have the prescription 



(27r)3 Vol 53 

where di represents the scalar degeneracy df = (/ + 1)^, or fermion degeneracy , df = 
Z(/ + 1). Thus, the one- loop corrected effective potential is = V^i + Vcw + K,hermai where 
V^i is the classical potential and 

11 °° 

i 1=0 
1 1 °° 

i 1=0 

The sum over / in the thermal part thermal is, as usual, convergent. However, the sum 
over / in Vcw diverges. Following [6], we regularize it by putting a momentum cutoff in 
the maximum allowed energy: 



/=o 1=0 

where 9{x) is 1 for x < 1 and for x > 1, and ui is the energy of the (scalar or fermion) 
mode. Note that the cutoff is on the energy, not on /, so the 9 function is slightly different 
for scalars and fermions (but such that that the energy of the last mode to be included 
in the sum is the same for both scalars and fermions). Then we replace the sum over / by 
integrals by using the Abel-Plana formula 



E = I 'dx Fix) + ^FiO) -2 Tdx ^^^^^ . (A.12) 
1=0 ^ Jo e -i 



The divergence in the sum over / is reflected in the first term. The integral will give 
rise to a A-^R^ piece which cancels between scalars and fermions, because this term is 
simply multiplied by the number of degrees of freedom which is the same for scalars and 
fermions. Then there is a A^R term. In fiat space, this term cancels out due to the 
vanishing supertrace formula. On (as noticed in [13]) it does not cancel out. However, 
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its coefficient is a constant independent of /x and independent of v. Therefore it will not 
affect the vacuum dynamics (which depends on the difference of free energies). Finally, 
there is the expected term proportional to log A. In a renormalizable theory, this can 
be canceled by the addition of a suitable counterterm. In the present case, we have an 
effective field theory model which is valid up to the energy scale A set by the dimensionful 
parameter (analogous to Fermi energy in the Fermi theory of weak interactions). 



B Kinetic term matrices 

We write the quadratic terms in the fluctuation Lagrangian as $''"05$ for the scalar part, 
and as "^^Op"^ for the fermion part, where 

$t = (0;(p), 0;(p), </,,.(-p), ^^(-p)) , (B.i) 
= ii^Up),i^Up),i^liip),i^Up),i^xi{-p),^x2{-p),^yi{-p),^ . (B.2) 

Consider first the Lagrangians (13.81) . (I3.10p . describing the model where the chemical 
potential is coupled to a U{1)b current. When = 0, Vy = 0, the kinetic term matrices, 
Os and Op, take the following simple form: 

Os = ^ diag [ {{u + iiif + + + Ag^^l) , ((w - iiif + + + %^A^) , 

((a; - l^^f + + + Ag^Al) , {{u + i^if + + + Ag^^l) ] (B.3) 
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We omit the (long) general expressions with v^^Vy 7^ as these are not used in the 
discussion. 

For the Lagrangian (I3.29P and (I3.30p . corresponding to a chemical potential coupled 
to a f/(l)R symmetry, we get the following expressions 
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Here {A, v} stands for {A^., v^} or {A^, Vy} as applied to the kinetic matrix involving the 
X or the Y superfield components, as the two chiral superfields X and Y are decoupled. 
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